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Dynamics and morphology of hole growth in a film of power hardening viscoplastic solid (yield 
stress ~ [strain-rate]™) is investigated. At short-times the growth is exponential and depends on 
the initial hole size. At long-times, for n > |, the growth is exponential with a different exponent. 
However, for n < |, the hole growth slows; the hole radius approaches an asymptotic value as 
t — > co. The rim shape is highly asymmetric, the height of which has a power law dependence 
on the hole radius (exponent close to unity for 0.25 < n < 0.4). The above results explain recent 
intriguing experiments of Reiter, Phys. Rev. Lett., 87, f86f0f (200f). 
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Dewetting of largely glassy polymer films below or near 
their glass transition temperature is a highly intriguing 
phenomenon in that the usual mechanisms of shape 
change such as viscous fluid flow, surface diffusion, evap- 
oration/condensation arc not dominant |Q . Thus, near 
glass transition temperature, solid-like plastic yielding 
and resulting flow of the glassy polymer appears to be 
the only mechanism of dissipation during dewetting. 

Indeed, Reiter recently observed the growth of holes 
in glassy polymer films with features very unlike those 
predicted/observed in purely viscous |, | films. The key 
features of the experiments are (i) the radius of the hole 
increases with time in an exponential manner (short time 
behaviour) (ii) the growth slows down considerably at 
long times, to the extent even of nearly stopping[Q (iii) 
the hole grows with a highly unsymmetric raised rim (iv) 
the height of the rim is nearly linearly dependent on the 
radius of the hole. 

In this Letter, we investigate the role of plastic yield- 
ing and resulting flow on the dynamics and morphology of 
growing holes in glassy polymer films. Our model is based 
on two physical premises: (i) the dewetting hole growth 
is driven by capillary forces, and (ii) the polymeric film, 
near its glass transition temperature, behaves as a strain- 
rate hardening viscoplastic solid^, i. c., the yield stress ~ 
(effective strain-rate)™, where n is the strain-rate harden- 
ing exponent (henceforth called simply as "hardening ex- 
ponent" ). The analysis uncovers a rich variety of regimes 
of hole growth some of which correspond to the experi- 
ments of Reiter |jj . 

The model consists of a flat film of initial thickness h 
with a hole of radius R a (see FIG. |l|); we focus on the 
growth of an existing hole rather than the mechanism of 
initial hole formation. As the hole grows, the radius in- 
creases to r a and the thickness becomes inhomogeneous. 
The aim of the analysis is to obtain the rate of the hole 
growth and the profile of the film including the height of 
the rim. 

Capillary forces are taken to be responsible for hole 
growth. As the hole grows there is a net reduction in 
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FIG. 1: Plane stress model of hole growth. The dashed con- 
figuration is the initial configuration of the film while the solid 
lines represent the configuration at time t. 



interfacial energy, the rate of which is given by 
E s = 2ir \S\ r r 



(1) 



with \S\ = \j S a — 7s/ — 7/a| is the spreading coefficient, 
where j s f is the interfacial energy between the film and 
the substrate, 7 sa is the surface energy of the substrate 
in contact with air and 7/ a is the surface energy of the 
film in contact with air. As the film grows, the height 
h(r Q ) of the rim increases and so does the surface area of 
the cylindrical surface of the hole. The rate of increase 
in the surface energy is given by 



E f = 2ir-i fa {h(r ) + h'(r )r )r 



(2) 



(We neglect the rate of increase of surface area on the 
top surface of the film, and this slightly overestimates the 
hole growth velocity). Applying principle of conservation 
of energy, we get 



E s — Ef + D p 



(3) 



where D p is the rate of plastic dissipation in the deform- 
ing film. 

The plastic dissipation D p (see final expression ( |l4| ) 
below) is obtained as follows. It is assumed that the film 
can slip freely on the substrate, i. e., the shear stress at 
the film substrate interface vanishes - this assumption 
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is justified by the fact that the fihn|l]] is deposited on a 
monolayer of PDMS. Neglecting inertial effects, dynam- 
ics provides a relation between the stress components (a 
plane stress assumption with a vanishing component of 
stress normal to the substrate (a zz = 0) is made) 



da rr 
dr 



= 0. 



(4) 



where a rr and agg are the radial and circumstantial com- 
ponents of the Cauchy stress tensor. Nonlinear kinemat- 
ics of large deformation is treated; the radial coordinate 
of the undeformed film is described by R and that of the 
deformed film by r. The deformation map can be ex- 
pressed as r = r(R, t) i. e., a point which was at a radial 
distance of R at the initial instance will move to a radial 
position r at time t. In particular, r a (t) = r(R ,t) is the 
radius of the deformed hole. The radial velocity is the 
time derivative of r defined asu = dr(R, t)/dt = r(R, t); 
the radial velocity of expansion of the hole is denoted as 
f = r(R ,t). The stretching rates e rr , egg and e zz are 
expressed in terms of the velocity field as 



dv v 
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-(e rr + egg). 



(5) 



Evidently, the kinematics ensures that the resulting de- 
formation is incompressible. 

We neglect elastic strains and treat the deformation as 
completely plastic. The polymeric material is taken to 
be a rate-sensitive plastic material (viscoplastic material) 
with power law strain rate hardening. The effective stress 
and the effective plastic strain rate at a material point 
are related by 



ke r ' 
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On substituting the above information into the equi- 
librium equation (H), we get, 



dr 
where 
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dv v 
dr r 



(10) 



The solution of (^J) provides an expression for the radial 
velocity 



v(r) = r = t I — 
r n 



(11) 



This equation can be used to obtain an expression for the 
thickness profile of the film 
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To calculate the plastic dissipation, we substitute all 
the relevant quantities into the relation 



ke n+1 h(r)rdr 



to obtain 



D p = 2nakh r 2 ( ^) F n (^- 



n+l 



(13) 
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with 



where a e 



2 Si o Si o with Sij 



&ij - \&kkbi 3 is the 
deviatoric stress tensor, e is the effective plastic strain 
rate defined as e = ^J^^ij^ij, k is a material constant and 
n is the hardening exponent. The polymeric material is 
assumed to obey an associative flow rule[|), i. e., 



e~ij — \Sij 



(7) 



where A is a number; the associative rule ensures that 
the plastic strain rate is normal to the yield surface. Re- 
lations (0), (|J), (H|) and (0) provide a complete set of 
equations for investigating hole growth in the polymeric 
film. 

Expression (0) can be used to express the nonvanishing 
stresses in terms of the stretching rates as 



2e rr + egg 
A : 
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(8) 



to find j = 



These expressions can be substituted in 
|fce n_1 - the stresses are thus expressed completely in 
terms of the stretching rates. 
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F n (x) = 2 Fi 



l-n l+3n~ . 2+3n+3rT . 
l+3n ' l+3n ' l+3n ' 



where 2-F1 is the hypergeometric function. 

Using (g), @ and @ in (§), the rate of hole growth 
is obtained as 
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(15) 



where t is an intrinsic time constant of the 
film/substrate system defined as 



to = 



and 








2n \S\ R ' 



(16) 



(17) 
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The rate of hole growth is controlled by two parameters 
(3 which is a measure of capillary forces, and n which is 
the strain-rate hardening exponent. 

The solution to ( |l5| ) can be investigated analytically 
in two limits, t < f and t t a . In the short-time limit 
(t -C t ), (p~5|) reduces to 



Ro to Ro 



(18) 



since F n (\) = 1; thus the hole growth is exponential for 
short times. It is clear that for the hole to grow the 
parameter (3 must be smaller than unity which provides 
a condition on the initial radius of the hole that allows 
for hole growth, i. e., (3 < 1 implies 



Ro > Rc 
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(19) 



where R c is the critical radius for hole growth ((3 = 
Rc/Ro)- If the initial radius of the hole is less than R c , 
the hole does not grow. As in a fluid film, the analy- 
sis for critical hole radius is most relevant when the film 
is sufficiently thick to ignore the excess intermolecular 
interactions. 

The long-time behaviour (t 3> t Q ) of the film is inves- 
tigated in two different regimes of n. First, when n > |, 
the rate of growth (TH]) for t ^> t Q reduces to 



To 

Ro 



Ro 



(20) 



(L n = lim-r^oo F n (x)~~) i. e., the growth continues to 
be exponential but with a different exponent L n /t - A 
plot of L n as a function of n is shown in FIG. |. Thus, al- 
though the short time exponent depends on the initial ra- 
dius R of the hole (through /?), the long time behaviour 
is independent of the initial size of the hole. Moreover, if 
the critical radius is small (R c <C h ), then the long time 
growth rate of small unstable holes ((3 <C 1) slows down 
compared to the initial growth rate L n < (1—0). 
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FIG. 2: Plot of L n vs. n. 



When n < i , the long-time behaviour is very different 
from the above case. Indeed, the hole radius does not 
become indefinitely large, but attains a asymptotic value 
r m given by 



(21) 



A straightforward analysis shows that the approach of r a 
towards r m is governed by 



1 -3n 



Ar a ^ 1 
Tm to \2ti F n (v m I Ro) 



where Ar Q = r m — r Q . Thus, r a approaches r m at a power 
law rate governed by the reciprocal of the hardening ex- 
ponent, 1/n. 
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FIG. 3: Hole radius as a function of time for various combi- 
nations of n and j3. 

The complete solution for ( |l5| ) for various combina- 
tions of n and (3 are shown in FIG. ||. All the features 
noted in the asymptotic analysis are reproduced. It is 
clear that for a given value of (3 the rate of growth is 
larger when n is larger, and for a given n, the rate of 
growth is larger when (3 is smaller. 

The time evolution of the film is also of interest; for 
r 3> r a , an asymptotic expression for the increase in the 
thickness of the film can be obtained (from (|l2])) as 



Ah(r, t) 



(1-n) (r (t)- 



Ro 2n ) 



(1 + 3n) 



(r»r D ) (23) 



FIG. |i| shows snapshots of film profiles for two different 
values of n for various times. A 3D view of the film profile 
is included in FIG. [5] for visualistion. 

Another quantity, important from the experimental 
point of view, is the dependence of the height of the rim 
on its radius. From (p"2"|), it is evident that 
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(24) 
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FIG. 4: Time evolution of film profiles. The profiles shown 
are for t/t a = 0.0, 0.5, f .0, f .5, 2.0. 




FIG. 5: 3D depiction of film profile, n = 0.25, (5 = 0.25, t/t a = 
5.0. 



which is, in general, a power law. However, if n is close to 
1/3, the exponent is approximately unity. In fact, when 
n = 0.25, h(r Q ) ~ r^ 5 and if n = 0.4, then h(r a ) ~ r° 75 
which implies that the relationship between the height of 
the rim and the radius of the hole is nearly linear in this 
regime. 

The main conclusions of this Letter are summarised as 
follows: (i) There is a critical value for the initial size of 



the hole. If the size of the hole is larger than this value, 
the hole grows, (ii) The short time growth velocity of 
the hole is exponential and depends on the initial size of 
the hole, (iii) The long time growth velocity of the hole 
depends on the strain-rate hardening exponent n. When 
n > i, the hole radius grows exponentially with time 
with a rate of growth that does not depend on the initial 
size of the hole. For n < |, the hole growth slows down 
and the radius of the hole attains an asymptotic value 
r m that depends on the initial size of the hole, (iv) The 
height of the rim is approximately linearly related to the 
radius of the hole when 0.25 < n < 0.4. 

The recent observations of Reiter [[l| on hole growth in 
polystyrene films near glass transition temperature can 
be explained by the above theoretical results, support- 
ing the premise that the hole growth in solid polymer 
films near T g is governed largely by yielding and subse- 
quent plastic flow which provides the dominant dissipa- 
tive mechanism. The theory also predicts a rich variety 
of growth behaviours which may help in the design and 
interpretation of new controlled experiments. Our theory 
for plastically deforming solids together with other recent 
theories of hole growth in viscoelastic fluids^, @, 0, [b]] 
and shear thinning viscous fluids^)\ provide complimen- 
tary models of dewctting over a wide range of tempera- 
tures from below T g to well above T g . 
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